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The Gaussian wavefunctional approach is developed in thermofield dynamics. We manufacture 
thermal vacuum wavefunctional, its creation as well as annihilation operators, and accordingly 
thermo-particle excited states. For a (D + l)-dimensional scalar field system with an arbitrary 
potential whose Fourier representation exists in a sense of tempered distributions, we calculate 
the finite temperature Gaussian effective potential (FTGEP), one- and two-thermo-particle-state 
QQ ■ energies. The zero-temperature limit of each of them is just the corresponding result in quantum 

^\ ' field theory, and the FTGEP can lead to the same one of each of some concrete models as calculated 

0^ . by the imaginary time Green function. 

The Gaussian wave-functional approach ( GWFA ) Q has become an important non-perturbative tool of quantum 
field theory and condensed matter physics g|| since Stevenson's advocation in 1985 [Q. It is simple, feasible, effective, 
and can reveal useful non-perturbative information (at least qualitatively) . A result from the GWFA can give a helpful 
[ guidance or be taken as a good start for a further investigation ||. Moreover, this approach is often used for some 
novel investigations or new ideas |6| . So the GWFA was inevitably generalized to finite temperature field theory, by 
a direct calculation of partition function from the first principle or by calculating the real-time or imaginary-time 
Green functions 0. However, thermofield dynamics ||, the third formalism of finite temperature field theory, has its 
own advantages over the other two formalisms : the closed time path and the imaginary time formalisms H. It can 
jJ^T reply to some questions which are not possible within the framework of the closed time path formalism, such as the 
structure of the thermal vacuum, the nature of the Goldstone states , and so on ||. Also, Thermofield dynamics can 
describe very naturally the time development for quantities near equilibrium || . Hence it is worth while generalizing 
£ — . the GWFA to thermofield dynamics. In 1986, through considering the A</>^ 0+1 ^ theory, I. Roditi gave the quantum- 
mechanical finite temperature Gaussian effective potential ( FTGEP ) within the framework of thermofield dynamics 
p0| . Recently, A. Mishra and H. Mishra calculated the FTGEP of the A<^ 3+1 ^ theory by the combination of the 
zero-temperature Bogoliubov and the thermal Bogoliubov transformation techniques ( TDBTT ) pj . In this letter, 
■ we intend to develop Roditi's pioneering work and establish the GWFA of thermofield dynamics ( TDGWFA ) in 
I \ (D + 1) dimensions. 

We consider a relatively general model with the Lagrangian density 
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£=-£W^x-n<M, (i) 

where cf> x = 4>(x), and the potential V(4> x ) has a Fourier reprensentation in a sense of tempered distributions p2| . 
5_i ' Potentials of many models, such as various polynomial models, sine-Gordon and sinh-Gordon models, have this 
property. We suppose that the field system Eq.(l) is immerged into a thermal reservoir with a fixed temperature 
T = ( fcfc represents the Boltzmann constant), i.e., we shall consider canonical ensemble. ( As for the Grand 
canonical ensemble, the generalization is straight foreward.) After manufacturing the thermal vacuum wavefunctional 
and the thermal particle creation as well as annihilation operators, we shall calculate the FTGEP, one- and two- 
thermo-particle energies, and then conclude this letter by a brief discussion. 

By the way, thermofield dynamics and the GWFA in quantum field theory will be not introduced in this letter, and 
good exposition of them can be found in Refs. §,§ and refs. respectively. 

Let us begin with quantum field theory. In the fixed-time functional Schrodingcr picture, the canonical conjugate 
operater to <f> x is H^ = — with the commutator [i^n^] = iS(x — y). For the convenience of renormalizing the 
theory of the system Eq.(l) later, we shall use the following normal-ordered Hamiltonian density with respect to any 
normal-ordering mass Q [ fl3| 

= + \(y^ x ) 2 - l -h[Q 2 ] + \Q 2 h[Q 2 ] +N Q [V{<t> x )] , (2) 



where V represents the gradient operator in £>-dimensional x-space, p = \p\, and I„[Q 2 ] = / p^jjp jp £+Q?)n ■ As for 
Mq\V{4> x )], making the Fourier transformation ( at least in a sense of tempered distributions [Q ) and using the 
Bakcr-Hausdorff formula, one can have 
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/DO JQ 
^Ly^yn^^i.m (3) 
-oo V27T 

Similar to the GWFA in quantum field theory 0^], we take following Gaussian wavefunctional as the quantum 
vacuum ansatz of T = field system Eq.(l) 



\tf >= N f eXp{~- I (<p x - <f x )f x y{(l)y - <Py)} , 

Z Jx,y 



(4) 



where ip x and f xy are variational parameter functions. Nf is some normalization constant, and depends upon f xy . 
Due to the translational invariance of vacuum, it is required that f xy = f(x — y). Also, in view of the symmetric 
form of the vacuum, we assume that f xy = f yx . Besides, the inverse f~ y of f xy has to exist. It is easy to show that 
tp x is the expectation value of the field operator with respect to the above trial vacuum \cp >. The creation and the 
annihilation operators to the vacuum Eq.(4) can be constructed as ]j]J| 

MP) = ( 2(2tt)W(p) )1/2 / e ^' S[fm<t>x ~ Vx) + lUx] (5) 

and 

4(p) = ( 2(27r) 1 g/(p1 ) 1/a I *"VM** - v.) - ^] (6) 

respectively, where f(p) is the Fourier transform of f xy . Of course, one can check the relations [Af(p), A\(j?)] = 
Sip 1 -p) and Af(p)\<p >= 0. 

In order to calculate average value of any physical quantity, thermofield dynamics needs the tilde system, the copy 
of the T = field system. Corresponding to the above, the tilde system has the normal-ordered Hamiltonian density 

= ^ + \{^k) 2 - \h[Q 2 ] + - A Q 2 h[Q 2 ] +M q [V{4> x )] (7) 



\(p >=N f exp{-\ I (<f> x - ip x )f xv (<f> y - <fy)} , (8) 



and the quantum vacuum 

. : ■--= \ r< .I'll I — ^ 

where fl x = igj- conjugates to the tilde field operator (f> x with the commutator [(p x , IL y ] = —i5(x — y). The appearance 

of the additional minus in the commutation relation is due to the rule of the tilde operation ||[|. Note that it is 
not necessary to distinguish the ordinary real function and average values of the T = system with those of its tilde 
system. From the rule of the tilde operation, the annihilation and the creation operators to the tilde vacuum Eq.(8) 
read 

and 

A *f® = ( 2(2^/(g) )1/2 I e ~**V®fc ~ W"> + ^ ■ ^ 

Of course, the relations [Af(p),A\(0')] = 8{p l — p) and Af(p)\<p >= hold. Additionally, each operator of the tilde 
system commutes with that of the original system. 

Basing on the T = field system and its tilde partner, now we take the following wavefunctional as the thermal 
vacuum ansatz of the T ^ field system 



\ip, /3 >= exp{i / [{<p x - ip x )g xy (f3)n y - Il x g xy (/3)(4> y - <p y )]}\(p<p > , (11) 

J X ,y 

where \iptp > represents the direct product of \cp > and \tp >, i.e., 
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\(f<P>=N}eXp{-~ I [{<j) X - (Py)f X y((f> y ~ (fly) + (4> X ~ ip x )f X y(4>y - Ify)}} . (12) 

Analogous with f xy , g xy ((3) = g yx {f3) = g(x — y,/3), and it is also a temperature-dependent variational parameter 
function to be determined. When it acts on the thermal vacuum wavefunctional Eq.(ll), any operator of the original 
system works only for quantities of the original system itself. So does that of the tilde system. According to Ref. [ JL4j , 
one can have ( here the Baker-Hausdorff formula is invalid ) 



exp{i [(<fr x - ip x )g xy (f3)Il v - II 

xQxy (/3)&, -¥>„)]} 

Jx,y 

= exp{- th(g(p, f3))n x (<t> y - tfy)} 

J xpy 

■exp{ / ln(ch(g(p, (3)))[{4> x - tp x )U y + (4> y - f x )U x ]} 

J xpy 

■exp{ / th(g(p, (3))(<f> x - ip x )fl y } , 

J xpy 

where f = i f xy J d ^(2*)° — > an< ^ 9(PiP) is the Fourier transform of the function g xy (0) . Hence one can write 
Eq.(ll) as 

\<p,P> = N 2 f exp{-~ [ [(cb x ch(g(p,P)) - 4> x sh{g(p, /?)) - <p x (l ~ sh(g(p,0)))) 

* J xpy 

■f{p){<j> y ch{g{f>,(})) - 4> y sh(g(p,/3)) - Vy (l - sh(g(p, /?)))) 

+$ x ch(g(p,p)) - 4> x sh(g(p,P)) -<p x (l- th{g(p,/3)) + sh(g(f, P))th(g(p, /?)))) 

■f(p)(4> y ch(g(p,0)) - cb y sh(g(p,[3)) - p„(l - th(g(p,p)) + sh(g(p, P))th(g(p, /?))))]} . (13) 

This wavefunctional is just Gaussian for 4> or for <f>. 

For the thermal vacuum Eq.(ll) or (13), one can manufacture the following annihilation and creation operators: 

Af(p,f3) = ( 2{2 ^ Dm ) V2 l{f(p)[e~ lfS ch(g(p, 0))(</> x - <p x ) 

-e^ s S h(g(p, 0)){4> x - if x )] + ie-^ s ch(g(p, 0))IL x - ie^ s sh(g(p, (3))fl x } (14) 

and 

A ^ P) = { 2(2n) D m )1/2 J{m\e**<*w, mt* - ¥>«) 

-e-^ s sh(g(p, (3))(4>x - <Px)] ~ ie^ S ch(g(p, (3))n x + ie~^ x sh(g(p, (3))Il x } . (15) 

with the relations [Af(p, /3), Aj (p 1 , (3)] = S(p l — p) and Af(p, 0)\ip, (3 >= 0. The thermo-particle creation and an- 
nihilation operators of the tilde system can be obtained from Eqs.(14) and (15) according to the rule of the tilde 
operation. 

Letting the thermo-particle creation operator (p, (3) act on the thermal vacuum Eq.(ll) or (13), we can construct 
the one- and two-thermo-particle excited states ( at the frame of the mass center ) as 

\l>=A\(?,0)\<p,p> (16) 

and 

|2 >= J d D pS(p)A\(p,f3)A\(-p,f3)\^(3> , (17) 

respectively, where S(p) is the Fourier transform of the S-wave function of the two thermo-particles with £(— p) = £(p). 
Of course, one can also manufacture the multi-thermo-particle states and the excited states of the thermal tilde system. 
Here, we discuss them no more. 
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Next we are ready to calculate the FTGEP. From Ref. B , the entropy operator is 

S = J dp[A f {p)A){p)ln(ch\g{p,m ~ A\(p)A f (p)ln(sh 2 (g(p, (3)))} , (18) 
and hence, the Hclmholtz free energy operator reads 

F = J Af Q [H x }-±S. (19) 

( Actually ksS is the entropy). Using functional integration technique and the integral formula J °° ^j=e~ a2x2 dx = 1, 
we obtain the average value 

da c 2 T T , a 



<<p,(3\jV Q [V(<f> x )]\<p,j3>= e—rV{-y/JM - h(Q 2 ) + <p x ) (20) 

J-oo 2 V 7T 2 

with the notation J n (g) = J ^^(/(gjP"^ ' ^ nus ma kmg some functional integrations leads to the vacuum average 
value of the free energy 



F[<p, g) f;lS\ = <<p,l3\F\<p,/3> 

\Ma) + ^(V^) 2 + \J ^ D J ch{2g{p,f})) - \h{Q 2 ) + -Q'U Ut) 



da 

[ch 2 (g(p,mn(ch 2 (g(p,m-sh 2 (g(p,mn(sh 2 (g(p,m} , (21) 



e-^V{~^J 1 {g)-h{QZ) + V > x ) 



0J (2tt) 

which is a functional of the functions <p x , <?, f. When the field system Eq.(l) has a thermal equilibrium with the heat 
reservoir at T ^ 0, its free energy is absolutely minimized. So minimizing the free energy with respect to f(p) ( the 
Fourier transform of function f xy ) and taking (p x — const. — ip, we obtain 



m = f(p, p) = vp^+jfiM) , (22) 

where 

/CO J 2 

e-^rv^H -VJi(ff)-/i(Q 2 ) + <p) (23) 
- oo W~K 2 

with V^ n \z) = d — f_^ c -^=(iVl) n V (Vl)e inz . Thus, we see that the variational procedure has enforced the 
function / depend both upon <p and upon j3 through the function g. 

Furthermore, regardless of the relation Eq.(22) and independent of the last variational procedure, one can minimize 
the free energy Eq.(21) with respect to g(p, (3) and obtain 

9(p,f3)=9(p,/3) = ln( eiWp ^ i ). (24) 

( Here, in the minimized result, we have made f(p, (3) take the place of f(p)). Now the function g is evidently related 
to /?, owing to the variational procedure. This point is similar to Ref. 11 1, but different from Ref. JTcp , where the 
expression between 9 and (3 is artificially given ( 9 there corresponds to g here ) . Substituting Eqs.(22) and (24) into 
Eq.(21), one has the FTGEP of Eq.(l) 



Vt(m> & <P) = -^(Mj ft V) = FYP = const., g xy ((3), f xy ; 0] 



l -Ug) - ^-Mg) - \l,{Q 2 ) + \Q 2 h{Q 2 ) 
da 



,e-^V(-^J x {g)-h(Q*) + V ) 



-oo2V^ 2 

\ J J^jD^ h2 (9(P^))Hch 2 (g(p,f3))) - sh 2 (g(p,(3))ln( S h 2 (g(p,P)))} . (25) 
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In this equation, "— >" means that /i takes the place of fi(ip,(3). Because of the nature of the minimizing procedure, \x 
should be chosen from the non-zero root of Eq.(23) and two end points of the range < /i < oo so that Vt[a*i f] IS 
an absolute minimum. Besides, sometimes the non-zero solution of Eq.(23) is multi- valued, and so in that case, the 
suitable root should be decided according to the stability condition 

(cV 2 ) 2 

= g J 2 (ff)[l + g U9) J ^ ^=e-^V^(^J 1 (g)-h(Q^ + V )} > . (26) 
For terms comprising J n (g)'s and I n (Q )'s in Eqs.(23), (25) and (26), we can write them as 

\j {g) \^,h{g) \l {Q 2 ) + \Q 2 h(Q 2 ) 
= \[h(p 2 ) i (Q 2 )} \v?h{» 2 ) + \q 2 Ii{Q 2 ) + Co 



T° 2 > 



Ji(g) - h(Q 2 ) = hip 2 ) - h{Q 2 ) + 2Ci , 

and 



\ I J^i ch2 (9(P,mn(ch 2 (g(p,m-sh 2 ( g (p 7 mn( S h 2 (g(p 7 P))^ 



J (27T) 

where C n = J r|^D ^f(p%'-i ■ Obviously, it is perhaps hard for Cq,C\ and Ci to have simple analytic expressions, 
but they are all finite. Moreover, for the case of D < 3, although i„(/x 2 (/3, ip))'s and I n (Q 2 )'s are divergent, the 
divergences can all cancelled each other in the above three expressions. So when D < 3, Eqs.(23),(25) and (26) have 
not any divergence, and accordingly no further renormalization procedures are needed. This indicates that the finite 
temperature field theory with Eq.(l) is renormalizable for the case of D < 3, at least within the framework of the 
FTGWFA. 

Now we are in a position for calculating one- and two-thermo-particle energies. For thermo-particle excited states, 
thermofield dynamics uses the following Hamiltonian H of the combined system of the T = field system and its tilde 
partner 

A/q [H] = [ {Af Q [H x ] - A/q [H x ]} , (27) 

J X 

and each thermo-particle excited state is the eigenstate of Nq [H] . Similar to the calculation of the FTGEP, using 
functional integration technique, one can obtain the one-thermo-particle energy 

m l(y ,/? ) = <1| ^ 1> = v /^+^/3) (28) 

and the two-thermo-particle energy 

_ <2|A/q[g]|2> 
™ 2 = < 2|2 > 



where 



fd°p[Z(p,(3)}* 

/°° A 2 
x 2^ e ~* rVin){ 2^ Jl{9) ~ h{Q2) + Z) ■ (30) 
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Minimizing ni2 with respect to S(p, (3) leads to a second kind of Fredholm integral equation about E(p, f3) Jig] , and 
solving it, we have 



£(p,/?) 



f dpch 2 (2g(p,/3)) 

C l ch(2g(p,(3)) ■ J 



/(p,/3)(2/(p, J 8)-m 2 ) 



, ^ 4 (y,ffl f dpch(2g(p,p)) 
1 + 16(2tt)° J p(p,/3)(2f(p,(3)-m 2 ) 



(31) 



with C some normalization constant. Substituting the last expression into Eq.(29) and noting the normalization of 
C, one can compute the energy 7712. From Eq.(28), one can see that [i(tp, (3) is just the mass of one thermal particle. 
Thus, Eq.(23) can be used for calculating the mass of a single thermo-particle. Moreover, the two terms in Eq.(29) can 
be interpreted as the kinetic energy of the two constituent thermal particles and their interacting energy, respectively. 
One can perhaps dicuss thermo-particle bound states and scattering states with the help of Eq.(29) |3| (1993). 

In this letter, we have developed the GWFA of thermofield dynamics in (D + 1) dimensions, and calculated the 
FTGEP, one- and two-thermo-particle-state energies of the system Eq.(l). Equations (25), (28) and (29) with some 
relevant relations are convenient for an investigation of any concrete model involved into the class of Eq.(l) ( including 



some bosonized models in condensed matter physics |16 ). When Eq.(25) is applied to A0 4 and </> 6 models, one can 
obtain the same, corresponding results as in Ref. ( Roditi, Okopinska, Hajj and Cea ). Again, contrasting the 
results in Ref. JTtJ with ones here, and noting the above-mentioned analysis about J„(g)'s and J n (<5 2 )'s, one can 
have the replacement rule between the FTGEP and GEP which indicated in Ref. ( J n {g) —> I n ([i 2 (<p))) . Still, the 
T = limit of each of Eqs.(25), (28) and (29) is identical with that in Ref. pj. Moreover, employing Eqs.(5), (6), 
(9) and (10), one can find that the exponential factor in Eq.(ll) is just field-operator form of a thermal Bogoliubov 
transformation, and therefore conclude that the TDGWFA is equivalent to the TDBTT. In fact, when I n (Q 2 )'s are 
deleted, Eq.(25) can give the non-renormalized equation Eq.(33) of Ref. jy]. Besides, if a non- uniform background 
field is considered, the term \{^f x ) 2 will be added to Eq.(25) with ip x replacing ip, and then one can have the 
following equation by T^l^&^sH = q 

^e-^Wi-y/Mg) - h{Q 2 ) + <p x ) = , (32) 
-oo 2V 7r 2 

the T = limit of which is identical with that of Ref. p7| . This equation can be used for considering quantum 
thermal solitons. Basing on thermal soliton solution, Eqs.(28) and (29) can perhaps give excited-state results of 
thermo-solitons. Of course, when the non-uniform background field is considered, the discussion about the difference 
between the BTT and GWFA in Ref. ji~7| is valid at the finite temperature case. Finally, the GWFA of thermofield 
dynamics can be generalized to the time-dependent case so as to investigate dynamics of statistical mechanical systems. 
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